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PACS 74 . 70 . Xa - Pnictides and chalcogenides 

PACS 74 . 20 . Rp - Pairing symmetries (other than s-wave) 

PACS 74 . 25 . Jb - Electronic structure 

Abstract -The pairing symmetry is one of the major issues in the study of iron-based supercon- 
ductors. We adopt a ten-orbital model by using the maximally localized Wannier functions based 
on the first-principles band structure calculations combining with the J1-J2 model for KFe2As2, 
the phase diagram of pairing symmetries is constructed. We find that the pairing symmetry for 
KFe2As2 is a nodal s^2y2 + s^2^y2-wa.-ve in the folded Brillouin zone with two iron atoms per unit 
cell. This pairing symmetry can explain the experiments observed nodes, and it also can be tested 
by future experiments. 



' , . The discovery of the superconducting iron-pnictide ma- 

^ ferial (1111-type i?eOFeAs, i?e=rare earth) [I] in 2008 

Qhas triggered great research interests which led to syn- 

Othesize similar iron-based superconductors, such as 122- 

' 'type BFe2As2 (B=Ba, Sr, or Ca) 0, Ill-type ^FeAs 

■^ ^A=alkali metal) [3], 11-tpye a-FeSe(Te) [1], and new- type 
J>Ka;Fe2-ySe2 [5] materials. Their crystal structures are 
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tetrahedral with the divalent iron square planes. A char- 
acteristic feature of these superconducting systems is that 
the band structure near the Fermi energy is derived from 
Fe-Sd orbitals with only modest hybridization of ligand-p 
drbitals [6]. Most of these compounds were reported to 
show superconductivity after doping or under high pres- 
sure and have the same robust tetrahedral layer structure. 

' To uncover the mechanism of superconductivity in these 
materials, the determination of pairing symmetry of the 
Superconducting order parameter is a good starting point. 
In most of the electron-doped and weakly hole-doped 1111 
and 122 compounds, the band structure calculations show 
that there are disconnected quasi-two dimensional (2D) 
hole and electron pockets [TlIS]. The former ones are cen- 
tered at the r point, while the latter ones are located at 
the M point of the Brillouin zone (BZ) of two iron atoms 
per unit cell. Strong scattering between the electron and 
hole pockets corresponding to a nesting wavevector [3|: 
q ~ (tt, tt), which could lead to the superconductivity with 
the so-called S^ pairing symmetry [5]. The order param- 
eter with this pairing symmetry is nodeless, and has a 



sign changed from the hole to the electron pocket. Such 
a S pairing symmetry is gaining increasing experimental 
supports [inilll] • For a heavily hole-doped compound like 
KFe2As2, the angle-resolved photoemission spectroscopy 
(ARPES) measurements [T21II3] indicated that the elec- 
tron pockets near M points were replaced by the ellipse- 
like hole pockets while the hole pockets near the F point 
became larger and closed to 7r/2 in the folded BZ of two 
iron atoms per unit cell. The superconducting transition 
temperature of this compound is Tc '^ 3.6 K which is 
much lower than those in other iron-based superconduc- 
tors. Under this situation q ^ {tt, tt) is no longer the nest- 
ing wavevector and the superconductivity pairing interac- 
tions could come from small q scattering. Although both 
experimental observations P^^ET] and theoretical investi- 
gations PM25] give nodal gap structures in KFe2As2, the 
pairing symmetry of this compound is still controversial. 
Specifically, experiments on small-angle neutron scatter- 
ing [131 indicates nodal lines perpendicular to the c-axis 
while experiments on thermal conductivity P^UTB] and su- 
perfluid density [TT] show a d-wave like nodal structure. 
Theoretically, based on the weak coupling approaches, a 
functional renormalization group study shown the dom- 
inant d-wave state [22] is most favorable, while Fermi 
surface (FS) restricted random phase approximation type 
spin-fluctuations analysis |23| and an analytical study in 
which the interactions are approximated by their lowest 
angular harmonics jMlIll] have revealed that s-wave and 
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Fig. 1: (Color online) (a) The schematic crystal structure of 
KFe2As2; (b) Sketch of the lattice structure of FeAs layers. The 
As ions • and • are located just above and below the center of 
each face of the Fe square lattice, respectively. 



d-wave pairing amplitudes have near-equal strength. An 
s-wave state vifith accidental line nodes that run either 
parallel to the c axis is also proposed by S. Maiti et al. 

m- 

In this paper, we address the pairing symmetry in 
KFe2As2 from a strong coupling approach, and treat the 
competing pairing channels on the same footing. Based on 
this approach, it has been shown that the pairing symme- 
try is determined mainly by the magnetic exchange cou- 
pling, as well as the FS topology for other iron-based ma- 
terials P7^l30j . In this sense, we construct a ten-orbital 
model by using the maximally localized Wannier func- 
tions based on the first-principles band structure calcu- 
lations and combine with the J1-J2 model to discuss all 
possible pairing symmetry. The mean-field phase diagram 
of a t-Ji-J2 model with correct band structure for this 
compound is obtained by performing the self-consistently 
calculations. We find that only the s^2^y2 +Sj.2j,2-wave can 
give rise to the nodal like structure in the folded BZ with 
two iron atoms per unit cell, and only this wave pairing 
symmetry can explain the experiments which evidenced a 
nodal order parameter [TSHIZ]- This result can be tested 
by future experiments, such as APRES measurements. 

In the first-principles calculations, the plane wave ba- 
sis method is implemented in the vasp code j31j . and 
the Perdew-Burke-Ernzerhof exchange correlation poten- 
tial ini] has been used. A 500 eV cutoff in the plane wave 
expansion and a 12 x 12 x 10 Monkhorst-Pack k-grid arc 
chosen to ensure the calculation with an accuracy of 10~^ 
eV, and all structures (lattice constants as well as inter- 
nal coordinates) were optimized until forces on individual 
atoms were smaller than 0.005 eV/A to obtain sufficient 
accuracy throughout the calculations. The crystal struc- 
ture is shown in Fig. HJa). 

Firstly, we focus on the electronic band structure, as 
shown in Fig. [2] It includes the calculated band structure 
(a), the corresponding FS (b) and the 2D cross-section 
of the FS at kz=0 (c) in the nonmagnetic (NM) state in 
the body-centered tetragonal unit cell. These results are 
similar to the experimental observations of this material 



Fig. 2: (Color online) (a) Energy band structure in the NM 
state of KFe2As2, and (b) the corresponding FS sheets crossing 
the Fermi energy in the band structure; (c) 2D cross-sectional 
{kz=0) representation of FS of the band structure calculations. 
The Fermi energy is set to zero. 



P^ l fn i lTT p g p g ] . As we can see in Fig. H there are six FS 
pockets, in which the three hole pockets (two of them are 
degenerate) are very cylindrical, suggesting a strong 2D 
behavior. However, there is another pocket near Z point 
with highly three-dimensional characteristics but vanish- 
ing around T point. In addition, four small ellipse- like 
hole pockets located at around the X points of the BZ. 
This result is consistent with ARPES results and previous 
first-principles calculations [T2lfT3] . Although these results 
have been reported, we present them here to facilitate the 
discussion on the orbitals characters and the possible pair- 
ing symmetry in superconducting state based on this FS 
topology. 

In order to clarify the orbitals characters in each pocket, 
we project the bands onto the five Fe-Sd orbitals shown in 
Fig. [31 From Fig. [U we can see that there is a dou- 
ble degenerate cylindrical hole pocket at F point with a 
dominant FeSdxz and Fe-Sdyz orbitals characters. How- 
ever, the Fe-3dz2 and Fe-3dxy orbitals characters become 
stronger around the Z point, which is slightly different 
from other iron-based materials j34j . The inner cylin- 
drical hole pocket is mainly contributed by Fe-3dx2-y2 
around zone center. It is interesting to point that the 
three-dimensional pocket around Z point is mainly coming 
from the As-Apz orbital character (partial Fc-3dz2 orbital 
character), which is similar to the one in FeTe [35] and 
KFe2Se2 121] ■ The orbitals characters of band structures 
exhibit the typical characteristics of layered structures. 

We now turn to discuss the superconducting pairing 
symmetry in KFe2As2. Because of the the strong hole 
doping in KFe2As2, the electron pockets at M point are 
replaced small hole pockets, and the hole pockets at F 
point are very large in the folded BZ of two iron atoms 
per unit cell, as shown in Fig. [Ijb). The nesting between 
r point and M point is absent. Therefore, we start with 
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Fig. 3: (Color online) Band structure of KFeaAsa, decorated pjg 4. (QoI^j. inline) The paramagnetic density-functional 

with partial characters of the Cg [top, d,2_y2 (a) and d,2 (b)] theory band structure (red dot lines) and a Wannier fit (blue 

and t2g [bottom, d^y (c) and d^^ + dy^ (d)] of the Fe-3d bands, ijjjgg) ^f ^^e ten orbitals in the vicinity of the Fermi surface 

The Fermi energies are all set to zero. o^^^ ^^^ Pg_3^ orbitals. The Fermi energies are all set to zero. 



a i- Ji-J2-type model: 



atoms per unit cell, 



H 



^la'^jvcr 



J2 ^m"«a- +J2J2 *u"4 
+J2 Yl (-^v • ^ji^ - 4"v"j>) 



J 



(1) 



where i,j denote the sites in the square lattice of FeAs 
layers [see Fig. [HJb)] and fi, v the orbitals, and t^'' is 
the transfer energy obtained from the maximally localized 
Wannier orbitals p7H39j . As shown in Fig. SI the ten- 
orbital tight-binding band structure reproduces the den- 
sity functional theory band structure rather accurately. 
The last two terms represent exchange interactions be- 
tween Fe-3d electron spin with parameters J\ and J2 de- 
noting intralayer nearest-neighboring (n.n.) and next- 
nearest-neighboring (n.n.n.) sites exchange interaction, 
respectively. The symbols (i,j) and {{i,j)) denote the 
summation over the n.n. and n.n.n. sites, respectively. 



The 5,^ = i E 



a,a' '^ifj.a'-^c^'^' 



and Tii 



E, 



-i/idr"-*M'^ 



Al,^(fc) = -^Xl'^fc'^fc' < C_fe/^;Cfc'^t >' 



J-: 



A2,^(fc) == -^Y<t>l(t>i' < C-fc>iCfcVt > 



(3a) 



(3b) 



Here the n.n. 
factors should have the forms (j)i = 4 cos ^ cos ^ (sx^y^ 



wave) or 4 sin 



2(cosA:a; -I- 



with a representing the Pauli matrices operating on the 
spin indices. For the purpose to discuss the superconduct- 
ing pairing symmetry, we obtain the following BCS-type 
Hamiltonian within the mean-field approximation: 



H^ J2 £feA<-4p<xCfe.^+E^M(fc)(Ct^-fcM4 + ^-C-) (2) 

The superconducting order parameters satisfy the fol- 
lowing BCS equations in the folded BZ with two iron 



and n.n.n. intraorbital pairing symmetry 

1 A fc, ky 

k 

^sin^ (dj^j^-wave), and 0^ 
cos fey) (s3;2^j,2-wave) or 2(cosfca; — cosfcy) {d^^ -y2-wave) 
in folded BZ with two iron atoms per unit cell. Note 
that the dxy-wave and d2,2_y2-wave pairing symmetries 
are nodal and s^2y2-wave and Sj.2_|_y2-wave are nodeless 
for any small doping parameter. In a previous study of 
a multiorbital t-Ji-J2 model for iron pnictides [27H30] . it 
has been demonstrated that the dominant pairing symme- 
try is governed by the intraorbital exchange interactions, 
and the interorbital exchange interaction only introduces 
quantitative modifications of the phase diagram. There- 
fore to keep the analysis simple, we will only consider the 
intraorbital pairings. The pairing order parameter A^(fc) 
for the orbital /i is linear combinations of four intraorbital 
pairings Ai,^(fc). 

The above equations: Eq.([2]) and Eq.dS]) can be solved 
self-consistently, we obtain the quasiparticle dispersion 
spectra Ek^- The pairing gap matrix is determined by 
minimizing the ground state energy density |28[|29| : / = 
Ez.fc.M tI<^1 < c-k^.ih^n > P - EkA^kf. - e^), where 
e'^ is the eigenvalue of ten-orbital tight-binding model. 

When the effective kinetic energy term is absent, the 
Ji dominating Sa;2y2-wave and dxy-w&ve pairing symme- 
tries are indeed degenerate. Likewise, the S2.2_|_j,2-wave and 
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Fig. 5: (Color online) Zero temperature superconducting phase 
diagram in the i-Ji-J2 plane obtained by solving the effective 
model in Eqs. ([2]). 



(i2;2_y2-wave dominated by J2 are energetically degenerate, 
and all four pairing states mix together when J2 is com- 
parable with Ji. Taking the effective kinetic energy term 
into account, the degeneracy states will be lefted, and the 
phase diagram of pairing symmetries is shown in Fig. [51 
The phase on the left upper corner has a time reversal 
breaking phase with Sx2^y2 + idj;2_y2-waNe dominated by 
n.n.n. intraorbital pairing interaction J2. When the n.n. 
pairing interaction Ji increases and is comparable to the 
n.n.n. pairing interaction J2, the (ia;2_j,2-wave will be sub- 
stituted by S3.2j,2-wave and lead to a S2;2+j,2 -|- S3.2j,2-wave 
in the middle part of Fig. [5] When the n.n. pairing in- 
teraction Ji further increases, the S2,2y2-wave will be sup- 
pressed and substituted by djjy-wave and lead to another 
kind time-reversal symmetry breaking state Sx2y2+idxy- 
wave. On the left lower corner of the phase diagram, when 
n.n. pairing interaction Ji is dominating and has a very 
small value, the system favors a nodeless 33,2^,2 -wave state. 
From the view point of experiments, the nodes have 
been observed by serval groups [TIH^ . Therefore, we 
conclude that only the Sx2+y2 + s^^yi-wave can give rise 
to the nodal structure, and only this phase can explain the 
experiments observed nodes based on the aforementioned 
discussions. The order parameter in this mixed phase can 
be written as 



kr. 



A^(fc) = Ao[cos-^ 



ky 

■COS- 



- S{cos 



kx + COS ky)] 



(4) 



where the value Aq is a constant and 6 depends on the n.n. 
interorbital pairing interaction Ji and n.n.n. intraorbital 
pairing interaction J2. When the n.n. intraorbital pairing 
interaction Ji is comparable to the n.n.n. one J2, the 
gap zero points will develop. As shown in Fig. [51 when 5 
reaches to a certain value, the contour of the gap zeros will 



Fig. 6: (Color online) A schematic plot for the gap nodal struc- 
ture in the folded BZ with two irons per unit cell. The blue 
lines denote the FS topology of the ten-orbital model. The red 
line is the contour line for zero gap value of s^- 
order parameter in Eq.@ with 5 = 0.2. 



cross the FS topology, which leads to the nodal behavior. 
This result can be tested by future experiments, such as 
APRES measurements. 

In conclusion, we have constructed a ten-orbital model 
by using the maximally localized Wannier functions based 
on the first-principles band structure calculations and 
combined with the J1-J2 model, all possible pairing sym- 
metries are examined by solving the self-consistently BCS 
equations. The mean- field phase diagram of a t-Ji-J2 
model for this compound is constructed. We find that 
only the Sx2+y2 + Sx2y2-waMe pairing symmetry can give 
rise to the nodal structure, and only this wave pairing 
symmetry can explain that experiments which evidenced 
a nodal order parameter. This result can be tested by 
future experiments. 
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